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Abstract 

We have made simultaneous torsional oscillator and transverse NMR measure- 
ments (at ~ 165 kHz) on 3 He contained within aerogels with nominal densities 
of 1% and 2% of solid glass. The superfiuid transition is seen simultaneously by 
both techniques and occurs at a temperature which agrees semi-quantitatively with 
that expected for homogeneous isotropic pair-breaking scattering of 3 He atoms by 
strands of silica. Values obtained for the superfiuid density p s in the 2% sample are 
in reasonable agreement with those observed previously. Coupling of the torsional 
mode to a parasitic resonance prevented accurate determination of p s for the 1% 
aerogel. We have identified other resonances coupling to the torsional oscillations 
as sound modes within the helium/aerogel medium. 
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1 Introduction 



The discovery of superfluidity of 3 He in low density aerogel glass in torsional 
oscillator experiments at Cornell Universityfl] and its subsequent investigation 
in NMR experiments at Northwestern University [2] raised many interesting 
questions; for example: 

• Would the transition seen in the torsional oscillator be coincident with the 
onset of an NMR frequency shift? 

• How many superfiuid phases are there and what is the phase diagram? 

• What is the nature of the pairing? 
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• Is it possible to explain theoretically the effect of the aerogel on the prop- 
erties of the superfluid phase (s)? 

In an attempt to address these questions we decided to perform simultaneous 
torsional oscillator and NMR experiments on the same sample. 

In this paper we report the results of torsional oscillator measurements on 
aerogel samples with 1% and 2% of solid density. Analysis of the NMR data 
is still in progress so we do not report these results in detail here. 



2 Experimental method 

The aerogel was grown within a spherical glass container with outer and inner 
diameters of order 10 mm and 8 mm respectively as shown in Fig. 1. Two sam- 
ples have been studied with nominal densities of 1% and 2% of that of solid 
glass. The 3 He entered the cell via a glass stem with an internal diameter of 
about 1 mm. In the case of the 1% sample the spherical envelope was com- 
pletely filled with aerogel but for the 2% sample a small region of bulk liquid 
was present near the stem as indicated in Fig. 1. The glass stem was glued to a 
berylium copper capillary with stycast 2850 GT epoxy[3]; the capillary acted 
as the fill line and also provided the torsion rod for torsional oscillations of 
the cell. The torsional oscillations were excited and detected electrostatically 
using electrodes mounted on the low pass vibration filter shown in Fig. 1. The 
vibration frequencies for the 1% and 2% samples were of order 850 Hz and 
970 Hz respectively. A two-phase lockin amplifier was used to measure the 
output from the detection electrode; from the in-phase and quadrature signals 
the resonant frequency and bandwidth of the oscillator were determined. Com- 
puter control of the driving frequency and voltage ensured that the oscillator 
was driven close to its resonant frequency and at a user-specified amplitude [4] . 

We were also able to simultaneously make transverse NMR measurements on 
the 3 He at frequencies up to about 200 kHz; a full account of the NMR ex- 
periments will be given later. Temperature was measured using an LCMN 
thermometer mounted in a separate tower. The LCMN was calibrated against 
the transition temperature of bulk 3 He[5]; the superfluid transition of the 3 He 
in the thermometer was identified from the change in warming rate of the 
LCMN. As we were able to detect the bulk superfluid transition in both the 
thermometer and the glass sphere we were able to monitor and correct for the 
small temperature differences between thermometer and aerogel. 
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Fig. 1. Schematic diagram of the experimental setup. The large moment of inertia of 
the demountable Be-Cu cylinder between the upper and lower torsion rods together 
with the lower torsion rod provides a low pass filter preventing high frequency 
vibrational noise reaching the torsional oscillator. The drive and detection electrodes 
are situated underneath the Be-Cu cylinder in the regions shown, thus avoiding the 
need to attach electrodes to the glass sphere and also allowing the use of electrodes 
of larger area. 

3 Superfluid transition temperature 



Fig. 2 shows the superfluid transition temperature T ca of the 3 He in the aerogel 
at various pressures relative to the superfluid transition temperature T c of bulk 
3 He. The values of T ca are plotted as functions of the bulk coherence length 
which we define as 



For the 2% aerogel the values of T ca were obtained from the torsional oscillator 
frequency measurements as described in Sec. 4. For the 1% aerogel the values 
of T ca were those below which a shift in NMR frequency from the Larmor value 
was observed; as explained in Sec. 5. it was not possible to use the torsional 
oscillator frequency to identify the transition reliably at all pressures. 

To the extent that the aerogel behaves like a dilute impurity which subjects the 
3 He superfluid to isotropic pair breaking scattering the reduction in transition 
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Fig. 2. Superfhhd transition temperature for 3 He in aerogel at various pressures. 
The circles and triangles are for 1% and 2% aerogel respectively. The theoretical 
curves are fits of the data to Eq. (2) using values of the mean free path I of 3380 A 
and 2440 A for the 1% and 2% aerogel respectively. 

temperature should be given by the following implicit equation[6,7] 




where ip(x) is the digamma function and I is the mean free path associated 
with the scattering. We have adjusted / to obtain the best possible fits to the 
experimental data; the values required were 3380 A and 2440 A for the 1% and 
2% aerogel respectively. As has also been observed in other experiments[7,8], 
there is a small but systematic difference between the observed pressure de- 
pendence and that predicted by Eq. (2) although this could be evidence for a 
small pressure dependence of I. The values of I are close to those obtained in 
a naive calculation for collisions of 3 He atoms with silica strands of diameter 
3nm[7]. Such a calculation predicts that / should be inversely proportional 
to the aerogel density; that our fitted values of I differ by less than a factor 
of two probably indicates that the real densities of our aerogel samples differ 
somewhat from their nominal values. We note that other measurements of T ca 
for aerogels with nominal density 2% of solid show considerable variation. Our 
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Fig. 3. Resonant frequency and bandwidth of the torsional oscillator as a function 
of temperature at 10 bar pressure for 2% aerogel. 

values of T ca lie between those of Ref. [1] and those of Refs. [2] and [9]. 

One important deduction that can be made from the success of Eq. (2) in 
describing the depression of the superfluid transition is that this indicates 
that the Cooper pairing within the aerogel is of the same type as in bulk, 
namely p-wave. 



4 Measurements for 2% aerogel 

Fig. 3 shows measurements of bandwidth and resonant frequency of the tor- 
sional oscillator for 2% aerogel filled with 3 He at 10 bar pressure; the temper- 
atures, T c and T ca , of the superfluid transitions of the small region of bulk 
3 He and of the 3 He in the aerogel are indicated. Above T ca the temperature 
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Fig. 4. Bandwidth versus resonant frequency of the torsional oscillator for T > T c 
(filled circles) and for T > T ca (filled triangles). Extrapolation of the linear relation 
to zero bandwidth as shown enables the frequency fo(0) to be determined. Values 
of vq(T) for T < T C a. were obtained by using the measured bandwidth at T and the 
extrapolation of the T > T ca data. The inset (top right) shows values of i^o(0) as a 
function of p; the straight line fit has an intercept which agrees with the empty cell 
frequency within experimental error. 

dependence is associated entirely with the small region of bulk liquid; the 
aerogel and 3 He within it behave as though rigidly locked to the motion of 
the oscillator. The increase in resonant frequency below T ca signifies the de- 
coupling of the superfluid fraction within the aerogel from the motion of the 
oscillator. The superfluidity of the 3 He within the aerogel does not produce 
any additional dissipation; apart from the small peak in bandwidth just below 
T ca and scarcely visible increments in dissipation even closer to T ca , which are 
discussed further below, the dissipation is associated entirely with the bulk 
liquid region. There is no evidence that the oscillator frequency and band- 
width depend significantly on magnetic fields up to 51 Gauss, the value used 
in most of our NMR measurements. 

To determine the reduced superfluid density p s /p within the aerogel we use 



where v(T) is the measured resonant frequency at temperature T as indicated 



p s = u(T) - z/ (T) 
p u e - z/ (0) 



(3) 
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by the solid line on Fig. 3, vq(T) is the frequency indicated by the dashed line 
that would be expected at temperature T in the absence of superfluidity in 
the aerogel and v e the measured frequency of the cell prior to filling with 3 He; 
z/ (0) is the frequency to be expected for the case where the 3 He within the 
aerogel is rigidly coupled to the motion of the oscillator but the bulk 3 He is 
completely decoupled. The value of ^o(O) could be obtained from both the nor- 
mal state data at sufficiently high temperatures and the data just above T ca ; 



in both these cases the viscous penetration depth 5 = y2i]/p n u is small com- 
pared to the size of the bulk liquid region leading to linear relations between 
bandwidth and resonant frequency as shown in Fig. 4. Extrapolation of these 
relationships to zero bandwidth (= 5 — > 0) allows ^o(O) to be determined. At 
all pressures the values obtained from the two extrapolations agreed within 
an experimental uncertainty of about lOmHz. To obtain the value of vq(T) 
for T < T ca the extrapolation of the bandwidth vs frequency for the T > T ca 
data was used as shown by the dashed line on Fig. 4. Since the dissipation 
was unaffected by the superfluid transition within the aerogel the measured 
bandwidth at the temperature concerned allowed us to obtain a value for the 
'unshifted' resonant frequency. This method ignores possible departures from 
hydrodynamic behaviour within the bulk region at low temperatures; since 
the variation of vq(T) with temperature is very weak (see the dashed line on 
Fig. 3), this is unlikely to lead to significant error. The values obtained for 
vq(0) at different pressures are shown in the inset on Fig. 4 as functions of the 
density p of bulk 3 He at the pressure concerned; the straight line through the 
data corresponds to that expected for rigid torsional oscillations of a sphere 
of diameter 4 mm and density 0.98p. The extrapolation of the line at p = 
to a value close to the measured empty cell frequency confirms our belief that 
the aerogel in the 2% sample moves rigidly with the glass envelope. 

Values of p s /p at different pressures obtained using Eq. (3) are shown in 
Fig. 5. The behaviour is qualitatively similar to that observed by Porto and 
Parpia[l,8] also for aerogel of nominal 2% of solid density; the data from 
Ref. [8] at 29 bar are also shown on Fig. 5. The temperature dependence of 
p s / p is very different to that of bulk 3 He. It rises more slowly as T decreases 
through T ca and approaches a value substantially less than unity at the lowest 
temperatures; this latter feature is more apparent in the data of Ref. [8] than 
in our data. Our data share with that of Ref. [8] the character that close to 
T ca the values of p s /p at different pressures can be superimposed to a good 
approximation by a translation along the temperature axis. For T/T ca > 0.75, 
p s /p is well fitted by a temperature dependence 



An example of such a fit is shown in Fig. 6. The values of the exponent b at 




(4) 
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Fig. 5. Temperature dependence of p s /p for 2% aerogel at the pressures indicated. 
The data of Porto and Parpia at 29 bar are shown by the dashed line. 

different pressures are given in Fig. 7. Fits to the bulk superffuid density over 
the same range of reduced temperatures do not fit Eq. (4) well near T c (where 
the correct limiting behaviour is p s /p oc (1 — T/T c )) and produce significantly 
smaller values of b ~ 1.25 — 1.30. We note that Porto and Parpia obtained 
somewhat smaller values of b than ours from logarithmic fits (ln(p s /p) vs 
ln(l — T/T ca )); also their values of b increased with increasing pressure. We 
used our fits to Eq. (4) as an objective way of determining the values of T ca for 
the 2% aerogel; this proved to be a more accurate way than from the NMR 
frequency shift data. There is no evidence for a systematic difference between 
the transition temperatures observed by the two experimental methods as can 
be seen from Fig. 8 which compares NMR frequency shift data with the value 
of T ca obtained from the torsional oscillator. 

Finally in this section we discuss the peak in bandwidth which can be seen 
just below T ca on Fig. 3. This peak was seen at all pressures above 7 bar but 
was largest at high pressures where a corresponding feature in the resonant 
frequency produced the small glitches in p s /p on Fig. 5. The occurrence of the 
peak just below T ca suggests a mode crossing of the torsional oscillator fre- 
quency with an internal mode of oscillation of the cell which has a frequency 
increasing from zero as p s becomes finite within the aerogel. A second smaller 
peak in dissipation was observed even closer to T ca ; there was evidence also 
for small additional dissipation between this second peak and T ca . These ob- 
servations suggest the existence of several internal modes such as might be 
associated with sound waves in the aerogel. 
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T (mK) 

Fig. 6. Fit of p s /p to Eq. (4) at 19.7 bar. The inset shows the behaviour close to T ca ; 
the apparent value of a small negative p s /p at T > T ca indicates a small error in the 
value of vq(T) used in the determination of p s /p through Eq. (3). In both figures 
the fitted curve is shown with a dashed line. 
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Fig. 7. Values of the exponent b giving the temperature dependence of p s near T c . 
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Fig. 8. NMR frequency shift as a function of temperature. The value of T ca indicated 
was obtained from the torsional oscillator measurements. 

Following McKenna et al[10], we assume that the normal fluid is rigidly locked 
to the aerogel. The sound speeds, c, are then the roots of the equation 

c 4 - c\c\ + c\) + c\c\ + ^(c 2 - c 2 )(c 2 - cl) = 0, (5) 



where c a is the speed of sound in the aerogel in the absence of the 3 He and 
Ci, c 2 and c 4 are the speeds of first, second and fourth sound respectively. In 
3 He, C2 is very small and can be set to zero in Eq. (5). One of the roots of 
Eq. (5) goes to zero as p s — > and the sound modes in a sphere of radius 
a associated with this root are calculated in appendix A. The two lowest 
modes occur at angular frequencies 2.082c/ a and 3.342c/a and lead to mode 
crossings with the torsional oscillator frequency at values of p s /p given by the 
theoretical lines on Fig. 9. The fit with the experimental values was obtained 
by adjusting c a ; the value used is 174ms _1 . McKenna et al quote a value 
of c a ~ 100 ms -1 for aerogel of 5% of solid density. Our value seems higher 
than might be expected for 2% aerogel since c a is expected to decrease as 
the density of aerogel decreases[ll]; it is possible that there will be variation 
between different samples of the same nominal density and also that c a may 
vary with temperature. It is clear from inspection of Fig. 9 that there is a 
systematic difference in pressure dependence between the theoretical curves 
and experimental points. We have no explanation of this although we note that 
the experimental values at low pressures are not well determined because the 
peaks were very small and somewhat broader. We note that the existence of 
coupling between the torsional oscillations and the sound modes indicates that 
our experiment does not have perfect rotational symmetry about a vertical 
axis. 
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Fig. 9. Comparison between the experimental (circles and triangles) and theoretical 
(continuous and dashed line) values of p s /p at which mode crossing of the torsional 
oscillator frequency with sound modes within the sphere occur. 

We consider also the possibility that the intersecting mode is a Helmholtz res- 
onance with pressure oscillations within the glass sphere producing oscillatory 
flow through the fill line. If our spherical container were filled only with 3 He 
then the geometry of our fill line is such that the Helmholtz frequency would 
always be more than a factor of two less than that of the torsional oscillator 
even when the superfluid fraction is 100%. The introduction of aerogel will 
decrease the frequency and hence the intersecting mode is not likely to be a 
Helmholtz resonance. There is a possibility that the sound modes discussed 
above could be modified by flow through the fill line. In appendix A we dis- 
cuss a simple model which takes both the flow through the fill line and the 
spatial variation of pressure within the aerogel into account. The conclusion 
is that the fill line impedance is sufficiently high that flow through the fill line 
is unlikely to have had a significant effect on the sound mode frequencies. 



5 Measurements for 1% aerogel 



Fig. 10 shows the bandwidth and the resonant frequency as functions of tem- 
perature at 15.0 bar pressure for the 1% aerogel. The bulk superfluid transition 
is clearly seen in both the amplitude and bandwidth despite the absence of a 
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Fig. 10. Resonant frequency (continuous line) and bandwidth (dashed line) of the 
torsional oscillator at 15 bar pressure for the 1% aerogel. 

macroscopic bulk superfluid region inside the sphere. The superfluid transition 
in the aerogel as determined from the NMR frequency shift is indicated and 
corresponds to an upturn in the resonant frequency with decreasing tempera- 
ture just as for the 2% aerogel (Fig. 3). However, just below T ca , the behaviour 
of the torsional oscillator is dominated by coupling to another resonant mode. 
We do not believe that this parasitic mode is a sound mode like those observed 
for the 2% aerogel. Indeed it may not be associated uniquely with the onset 
of superfluidity of the 3 He in the aerogel since at low pressures the coupling 
to this mode appears to be evident above T ca ; it is this coupling above T ca 
which prevented identification of T ca from the torsional oscillator data at low 
pressures. The large magnitude of the coupling between the mode and the tor- 
sional oscillations together with the large bandwidth of the torsional oscillator 
at all temperatures in comparison with the 2% data suggest that significant 
shear motions of the 1% aerogel are being excited. Further evidence for this 
is the apparent discontinuity in the 'background' frequency of the torsional 
oscillator associated with the mode crossing. This possible discontinuity pre- 
vents the deduction of reliable values of superfluid density for the 1% aerogel 
at lower temperatures although we are currently searching for a theoretical 
description of the intersecting mode which might make this possible. 

The smaller intersecting resonance which can be seen at lower temperatures 
on Fig. 10 is we believe due to a sound mode like those observed for the 2% 
aerogel. If we ignore the apparent discontinuity in the background frequency of 
the torsional oscillator mentioned in the previous paragraph we can estimate 
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a value of p s at which this second resonance occurs and then by proceeding as 
for the 2% aerogel we can estimate the speed of sound c a in the 1% aerogel. 
The value obtained is 55 ms -1 which, in view of our neglect of a possible 
background frequency discontinuity must be regarded as an upper limit; the 
value is less than that for the 2% aerogel as would be expected. 

For the 1% aerogel we investigated the effect of the addition of small quantities 
of 4 He on the torsional oscillator. About 3% of 4 He was sufficient to replace the 
solid He layer on the aerogel surfaces as indicated by the absence of a Curie- 
Weiss contribution to the NMR signal strength[2]. The addition of 4 He caused 
a dramatic reduction in the size of the parasitic resonances and of the oscillator 
bandwidth at higher temperatures; coupling to the intersecting resonances was 
almost completely absent at high pressures. We do not have any explanation 
for this observation. Although the parasitic resonances were absent we were 
still unable to determine the superfluid density within the aerogel because 
the addition of the 4 He also introduced a large thermal boundary resistance 
between the LCMN and the 3 He which preventing us from measuring the 
temperature. 



6 Conclusions 

In this concluding section we return briefly to the questions posed in our intro- 
duction. Our experiments show that the superfluid transition of 3 He in aerogel 
as indicated by the appearance of a finite superfluid density coincides with the 
onset of a shift in the NMR frequency. Our torsional oscillator measurements 
provide no evidence for more than one superfluid phase and do not identify the 
nature of the pairing although our observation that the measured superfluid 
density is independent of applied magnetic fields up to about 50 Gauss might 
be interpreted as implying that the superfluid phase is one with an isotropic 
superfluid density. We are hoping that the completion of our analysis of our 
NMR measurements will provide further information on the questions of the 
number of phases and the nature of the pairing. The fact that we measure two 
different properties on the same sample of aerogel should provide a stringent 
test for theories purporting to explain the effect of aerogel on the properties 
of the superfluid. 
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A Sound Modes and Helmholtz Resonance 

We consider first the sound modes in a spherical cavity of radius a completely 
filled with aerogel containing 3 He. The small value of C2 for 3 He means that 
the effect of temperature gradients within the cavity can be ignored and the 
equations of motion which describe the 3 He/aerogel combination are then[10] 
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dp . . . 

— = -V • (p s v s + p n v„), (A.l) 

aV " = -ivp, (A.2) 



dt p 
<9v : 
~dt 
dp 



(Pa + Pn)^ = -^Vp-V Pa , (A.3) 

ar p 



01 -V-( Pa v n ), (A.4) 

where p and p a are the densities of helium and aerogel, p is the pressure 
acting on the helium and p a is the pressure acting on the aerogel due to its 
elastic distortion; the motion of the aerogel and normal fluid are assumed to 
be locked together. Assuming that small variations in p and p a are related 
to the corresponding densities by 5p = c\8p and 5p a _ = c^Sp^, we obtain 
from Eqs. (A.l) to (A.4) the following equation for small harmonic departures, 
5p = p'exp(icut), of p from equilibrium 

vt4 f , 2/1 ' Pn c l \ ^2 / i 4(Pa + Pn) / „ / \ r\ 

V p + u \ ^ + ^ + - 7r7 ^ 1 )V p + u p = 0, (A.5) 



where c\ = c\p s / p. 

The solutions of Eq. (A.5) appropriate to our spherical geometry are of the 
form 

p' = (Ajtihr) + B 3l (k 2 r)) Y lm (9, 0), (A.6) 



where A and B are constants of integration, ji(z) is a spherical Bessel function, 
Yi m {9, <p) is a spherical harmonic and we are using the notation of Abramowitz 
and Stegun[12]; the spherical Bessel functions yi{kir) and yi{k 2 r) can be ex- 
cluded from the solution because they diverge as r — > 0. The wave numbers 
ki and k 2 are the roots of 



^^(l + i + ^W%±£>=0, (A.7) 



and are thus the wavenumbers associated with the sound speeds which satisfy 
Eq. (5) in the limit c 2 — >■ 0. The boundary conditions to be applied to the 
solution are that 

v n • f = v s • f = at r = a. (A. 8) 



The radial components of the v n and v s can be found for solution (A.6) by us- 
ing Eqs. (A.l) to (A.4). Fortunately the boundary conditions produce no mix- 
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ing of the terms of different wavenumber and the modes occur at wavenumbers 
k (= ki or k 2 ) which satisfy 

= Jl _ l (ka)- l -±± Jl (ka) = 0. (A.9) 



The roots are thus characterised by two numbers n and /, where n specifies 
the number of radial nodes (excluding r = 0) in the variation of p. The values 
of ka (n,l) for the 7 lowest modes are 2.08157 (0,1), 3.34210 (0,2), 4.49341 
(1,0), 4.51411 (0,3), 5.64670 (0,4), 5.94037 (1,1), 6.75645 (0,5). The two lowest 
modes with n, I — 0, 1 and 0, 2 are the modes plotted in Fig. 9. 

We now consider the possibility that the modes may be affected by flow 
through the fill line. Since a spherical cavity on the end of cylindrical fill 
line is difficult to calculate we consider instead a simplified model: a spherical 
cavity with a fill line entering at the centre. We take the end of the fill line 
to be surrounded by a small spherical region of bulk liquid of radius b. Only 
modes of spherical symmetry will have a pressure variation at the centre of 
the sphere. For these modes the general solution of Eq. (A. 5) is 

p' = (Atfofar) + B m {k ir ) + A 2 j {k 2 r) + B 2 y (k 2 r)) . (A. 10) 



Because of the small region of bulk helium and the fill line at the centre 
of the sphere it is no longer possible to ignore the spherical Bessel function 
yo(z) = — cos(z)/z which diverges as z — > 0. The boundary conditions to be 
applied at r = a are Eqs. (A. 8) as before and at the boundary of the bulk 
liquid region, r = b, we require that the variations of p & should vanish and 
that the ratio of pressure variation bp = p' exp(iut) to mass flow out of the 
aerogel M = —4nb 2 (p s v s + p n v n ) should be appropriate to the geometry of our 
fill line 

sm = ! upL 

M Psbo- 



where L/a is the ratio of fill line length to cross sectional area and p sb is 
the superfluid density in the fill line. Eq. (A. 11) follows from the equation of 
motion 



for the helium within the fill line. Note that we ignore the variation of pressure 
inside the sphere of radius b. 
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Applying the boundary conditions to the solution given by Eq. (A. 10) leads to 
a rather complicated condition for determining the frequency of the spherically 
symmetric modes. The terms in Eq. (A. 10) of wavenumbers k\ and hi are no 
longer decoupled. We do not discuss the details here but report the general 
conclusion that the impedance of our fill line (L/a) is sufficiently high that 
the modes under consideration are essentially sound modes with little flow 
through the fill line. 
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